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^ ■ Abstract 

Tlie derivation of a new family of magnetic fields inducing exactly solvable spin 
evolutions is presented. The conditions for which these fields generate the evo- 
lution loops (dynamical processes for which any spin state evolves cyclically) are 
5r , studied. Their natural connection with geometric phases and the corresponding 

calculation is also elaborated. 

1. INTRODUCTION 

In nonrelativistic quantum mechanics, the pure states of spin-1/2 systems are usually 
represented by two-component ket vectors \ip) = {ipi,ip2), where '4'i,'4'2 £ C"*^ are the 
components of \ip) along the two orthogonal eigenstates |+) and |— ) of 5*2 with eigenval- 
ues h/2 and —h/2 respectively. If the spin is placed in a homogeneous time-dependent 
magnetic field B(t), then the Hamiltonian is described by H{t) = — /iB(t) ■ S, and the 
Schrodinger equation governing the evolution becomes explicitly time-dependent: 

^nf^\m) = Hit) m)) = -/iB(t) ■ S|^(t)). (1.1) 



In order to solve ( |1.1D , usually it is rewritten in terms of the time-evolution operator 
U{t), U{t = 0) = /. The most elementary case arises when [H{t),H{t')] = for all 
t 7^ t' [jl|. In the case when [H{t), H{t')] ^ 0, however, to sum up the continuous Baker- 
Campbell-Hausdorff exponent becomes hard |0], and thus it is not easy to construct 
exact solutions to (|1.1| ). Due to this, the approximate methods (see for example [H 
and references quoted therein) or direct numerical techniques 0] are the standard tools 
applied in the general case. Is there any optional technique leading to exactly solvable 
situations? The answer seems to come of the so called inverse techniques, which take 
full advantage of the geometrical picture involved in the spin-1/2 description 0. 

Suppose that somehow one knows the time evolution of the system, i.e., the state 
vector \ip{t)) satisfying (1.1) for all t, but there is no information about the external 
field B(t) driving the system. Can one find at least one field which could in principle be 
created and would dynamically induce that state? This is the essence of the dynamical 
manipulation problem whose main ideas were developed by Lamb and followed later 
by other authors 0, ^, Q. In the spin-1/2 case (1.1), this technique was successfully 
applied to generate exactly solvable situations 0. It would be interesting to show that 
this is so using the standard direct approach. 

The aim of this paper is to prove explicitly that given the family of magnetic fields 
derived in P, they lead to exact solutions to equation (1.1). Moreover, we shall show 
that those magnetic fields can induce some special dynamical processes on the system 
such that the evolution operator becomes the identity (modulo phase) at some t = t, 
i.e., U{t) = e^'^I. Such kind of processes have been widely studied under the name of 
evolution loops (EL). It has been also proposed that the EL can be perturbed in order to 
induce any unitary operator as a result of the precession of the distorted loop ||^ . This 
proposal has been successfully realized by inducing the squeezing inside of a modified 
Penning trap and by rigidly displacing the wavepacket in a magnetic chamber perturbed 
by homogeneous time-dependent electric fields 0. Here, the EL will mean just closed 
orbits on 5*^ (independent of the initial condition) which can be either periodic or 
aperiodic and whose corresponding geometric phases can be simply evaluated. 

The paper is organized as follows. In Section 2 we briefly sketch the derivation of 
the magnetic fields of [Q with a discussion of their basic properties. In section 3 we 



directly solve (|1.1| ) taking for B(t) the analytic expressions of these fields. Then we 
shall establish the conditions on the fields in order to induce the EL on the system. 
To conclude, we will show that these EL give place to integrable expressions for the 
corresponding geometric phases. 

2 . THE INVERSE METHOD 

For the set of pure states of spin-1/2 systems, a geometric picture is easily found by 
noticing that the inner product between any two ket vectors is invariant if both are 
multiplied by an arbitrary unimodular complex number A G C^, |A| = 1. Hence, the 
space of physical states corresponds to the projective space CP(2) which, as is well 
known, is usually modeled by the Riemman sphere S"^ with each point of the surface 



representing a pure state of the spin-1/2 system. 

Now, according to Eherenfest theorem the mean value of the operator S = ^a; i + 
Sy3+ Sz^, where Sk = (^/2) cr^ and [Sk, Se] = ihtkemSm, evolves as: 

|(S) = -/iB(t) X (S) = -b(t) X (S), (2.1) 

where b(t) = /xB(t). The identification n = (2/h){S), n ■ n = 1 shows that ( p.l| ) is 
precisely the dynamical rule governing the evolution on S"^. In the direct approach 
the initial vector n(0) and the field b(t) are given, and one looks for the solution 
n(t) to (2.1). Here we assume that the spin state n(t) is given and rewrite (2.1) as 
n(t) = M[n(t)]b(t). An element of arbitrariness arises by noticing that the matrix 
M[n(t)] is antisymmetric, and hence its determinant is equal to zero. Thus M~^[n(t)] 
does not exist, which does not allow to determine uniquely b(t). Let us take the third 
component &3(t) also as given; henceforth the other two components become: 

hit) = [b3{t)ni{t) + n2{t)]/Mt), b2{t) = [h{t)n2{t) - ni(t)]/n3(t). (2.2) 

Notice that, departing from a given point n(0) G S'^, any other point n(t) G S"^ can 
be achieved by a set of successive infinitesimal rotations encoded in R(t) G 5*0(3), i.e., 
n(t) = R{t)n{0). Therefore, the field b(t) given by (2.2) depends on the generic motion 
(a generalized rotation) and the initial condition. Thus, two paths with common R{t) 
but different n(0) are induced by two different fields with the same 63 (t). Would it be 
possible that trajectories sharing the same R{t) determine a unique b(t)? In order to 
find the answer, let us consider the case when n(t) rotates simultaneously around two 
fixed directions with variable angular velocities. By simplicity, one of these directions 
is fixed along k and the other one along a vector e^ on the x — z plane at an angle x 
from k, e^ = sin xi + cos xk. The rotation matrix is: 

R{t) = i?3(/5(t))i?x(«W) = i?3(/5(t))i?2~'(-x)^3(-«(t))i?2(-x), (2.3) 

where R2{'^) and Rslw) are finite rotations by 00 around j and k respectively and 
a(0) = /5(0) = 0. Using ( p.2|) , one will find a b(t) dependent of the initial condition; the 
field independent of n(0) arises after imposing the restriction b^it) + [5{t) = a(t) cosx, 
and then: 

b(t) =d(t)sinx[cos/?(t)i + sin/?(t)j] + [d(t) cosx - /?(t)]k. (2.4) 



Let us remark that, due to its dependence on two arbitrary functions, equation ( |2.4|) 
represents a wide family of analytically solvable fields inducing the rotation described 
by (2.3) on any initial vector n(0). An interesting result comes out by analyzing the 
following physical situation. Suppose that the spin initially points along an arbitrary 
direction n(0) = e+; what is the probability that at time t the spin will be in the 
corresponding orthogonal state e_? Notice that orthogonal vectors {ip) on the Hilbert 



space Ti correspond to antipodal points on the sphere S*^. By simphcity, let us choose 
e± = (0, 0, ±1); then the probabihty transition is given by: 

P+^-{t) = [1 - n^{t)]l2 = sin2xsin2[a(t)/2]. (2.5) 



As our treatment is exact, equation ( |2.5| ) is indeed a generahzation of Rabi's formula 
for any t and arbitrary a{t). It reduces to the standard Rabi expression when we take 
a{t) = aot for small t (see e.g., Rabi et.al. []TU| and Shirley |^). Therefore, by choosing 



specific forms for a{t) and /3{t) in (|2.4|) one is led to obtain diverse particular cases of 
b(t) 1^, some of which could have been previously discussed. 

Let us finish this section by remarking that the product bk(t)Sk of ( p..lD satisfy 
bk{t)bi{t') [Sk,Si\ = ihekembk(t)bi{t')Sm- Hence, for arbitrary functions a(t) and P(t) 
(see the expressions of bk{t) in ( |2.4| )), the Hamiltonians H{t) = — b(t) ■ S at different 
times do not commute, i.e., [H{t), H{t')] ^ 0, t ^ t' . Then, we have arrived at a family 
of exactly solvable Hamiltonians which, at first sight, should be solved by approximate 
or numerical methods when using the direct approach. 

3 . THE DIRECT METHOD 

In this section we are going to solve ( |1 . 1| ) with the magnetic field ( |2.4{ ) using ordi- 
nary quantum mechanical operator methods. The key point is that, by means of the 
commutation rules [ak, crp] = 2iekemO'm, and the expression 

g-*,3(t)<x3/2^^ e*/3(t)<x3/2 ^ cos/?(t) di + sin (3{t) as, 
the Hamiltonian H{t) = — b(t) • S can be rewritten as: 

H{t) = e-^^W-^^/^ [Heffit) + n${t)a3/2] e^^^'^^/^ (3.1) 

where Heffit) is defined by 

- {2/h) Heff{t) = d(t)(sinxcTi + cosx^a). (3.2) 

Now, let us introduce a new reference frame which rotates with angular velocity (3{t) 
around k. Notice that rotating frames are very useful in magnetic resonance because 
when successfully used, the original problem can be mapped into a static one, which 
simplifies the physical analysis of the problem |l^ . The key transformation in this case 
is given by U{t) = e"*'^^*)'^^/^ W{t), where W{t) is a new unitary operator satisfying the 
equation: 

in^^ = H^ff{t)W{t). (3.3) 

It is clear now that Heff{t) is the Hamiltonian in the rotating frame, and thus W{t) 
represents precisely the evolution operator in that frame. Notice that Heff{t) satisfies 



[Heff{t),Heff{t')] = 0; henceforth, the solution of ( p.3|) is given by just integrating 

Heffity. 

W(t) = e*-~*^^-'-^o ■'^^^'^'^*''''^*' = e*"*-*-'*'^''^^'^^^'^"'^^'^^^/^ (3.4) 

= cos(Q;(t)/2) + i(sinxo"i + cosxo"3) sin(a;(t)/2), 

where we have used again the algebraic properties of 0-^. Notice that at t = equa- 
tion (|3.4| ) reads W{0) = I. Hence, initially the evolution operators U{t) and W(t) 
coincide. At arbitrary times t 7^ 0, there will be a time-dependent factor operator 
exp(— ■i/9(t)cr3/2) making the difference between the descriptions at the lab and at the 
rotating frames. 

Let us construct now the generic expression for the ket vectors \ip{0)). First, the 
normalization condition gives |'?/'i(0)p+ |'?/'2(0)P = 1. It suggest the following represen- 
tation ipi{0) = cos(^o/2), ip2{0) = sin(^o/2). Now, the normalization remains invariant 
if we take '?/^i(0) = exp(— z0i) cos(6'o/2), ■?/'2(0) = exp(i02) sin(6'o/2). As mentioned at 
Section 2, the multiplication of \ip{0)) by a common phase factor A = e**^''^i~'^2)/2 (^oes 
not change the representative of the corresponding physical state. Hence, we can take: 

/ ^i(O) \ / cos(f ) e-^<^«/2 \ 

1^(0)) = = , (3.5) 

V ^2(0) / V sin(f)e^<^«/2 J 

with 00 = 01 + 02- We are using the half angle convention whose utility will be apparent 
below. The solutions \ip(t)) to ( |1.1| ), with the magnetic field b(t) given in ( |2.4| ), arise 
just as time displacements induced by the evolution operator U{t) acting on the initial 
vector (|3.5|) : 

where we have used (p.4|) . Let us remark that the map |'0(t)) -^ {2/h){S)(t) = n(t) 
reproduces the results derived in Section 2, i.e., the vector n(t), connected with \ip{t)) 
in ( p.6|) by this map, is the result of two simultaneous rotations performed by the vector 

n(0) = {2/h) (^(0) I S I V'(O)) = (sin Oq cos 0o, sin Oq sin 0o, cos ^0) (3.7) 

around the two directions and with the two angular velocities characteristic of the 
matrix (^.3[). Notice that, from a geometrical point of view, the above map corresponds 



to the Hopf map [0, [T^ , which is useful in the description of monopole magnetic charges 



in the fibre-bundle formulation of electrodynamics [0], and provides an interesting 



geometrical interpretation of the Aharonov-Bohm effect |14, O 



3.1. Evolution Loops and Geometric Phases 

Let us formulate now the requirements which has to be satisfied in order to induce the 
evolution loops. With this aim, let us rewrite the evolution operator in (|3.6|) as: 



f/(t) 



cos — -— — 20-3 sm — -— 



fa{t)\ ... , . fa{ty 

cos — -— + «(sm xo'i + cos X'^3) sm — -— 



It is apparent that a{T) = 2i7T and /3{t) = 2imi, with £, m G Z, produce U{t) = 
cos{m7i) cos(£7r) = cos(m + i)n = (—1)™+^/. The loop conditions are thus: 

a{T) = 2in, /5(r) = 2m7r, m,£eZ;r>0. (3.8) 

Notice that there is a strong loop condition in which the evolution operator becomes 
the identity sensu stricto. In our spin-1/2 case this strong loop condition consists of 
the restriction (3.8) and the additional requirement m + i = 2k, k E Z. Here and 
throughout the paper we will use the relaxed loop condition (3.8). 

The conditions (|3.8|) have been intuitively used in the case when the magnetic 
field (2.4) rotates with constant angular velocity, and then the vector n(t) describes a 
hypocycloid on S"^ [^ (similar results can be found in |jl6| and Zhang et. al. @]). That 



case is recovered here by taking a{t) = aot, (3{t) = (3ot and by forcing the loop condition 
(|3.8|) . Other selection of the functions a{t) and P(t) allows us to generate deformed 
versions of such a case and generalized nontrivial cases (see the specific examples of 

i). 

The loop condition (|3.8|) and equation (|3.6|) give \iP{t)) = (— 1)™'+^|'?/'(0)). As 



Aharonov and Anandan have shown, any cyclic quantum state has naturally associated 
a geometric phase 7 characterizing somehow the projective Hilbert space curvature (see 
also the general geometric treatment introduced long ago by Mielnik [|13). In the spin- 
1/2 case, it turns out that 7 = — Ar2/2, where Afl is the solid angle subtended by the 
oriented closed curve n(t): 

AQ= / ^^ ^—^dt. 3.9 



Let us remark that, although the general expression ( |3.9| ) has been recurrently 
studied and discussed in the literature [15-20], it is not always possible to perform the 
involved integrals (see [M and references quoted therein). Therefore, it is interesting 
to look for explicit expressions for 7. For the magnetic fields (2.4) we have gotten the 
generic time-evolution (|3.6| ), and this added to the loop conditions ( p. 81) allow one to 
simplify considerably the calculation: 



7 = [a — m + COS 



1 r ■ 
i(^o — x)("^cosx — £)]7r sinxsin(6'o — x) / P(t) cos a {t)dt, (3.10) 

2 JO 

where in ( p. 71) we have taken 0o = 0. Notice the arising of the atypical integral term 
in ( |3.10| ). This formula generalizes the corresponding expression for the traditional 
rotating magnetic field with constant amplitude and angular velocity [^ . 

In order to illustrate the generality of our expression ( |3.10| ) , let us consider some 
particular cases of (^^ ). The first simple case arises by taking 63 (t) = 60 in ( ^^) , and 
thus 

7 = [£ — m + cos(6'o — x)("^cosx — ^)]7r H sin x sin(6'o — x) / cos a{t)dt. (3.11) 

2 JO 



The remaining integral depends just on a{t) and r. Notice that it vanishes for the 
simplest time-dependent function a{t) = aot with the loop condition aoT = 2iTT, which 
is indeed the case discussed in [|15| and [16|. The first nontrivial case of ( p.ll ) 



arises 



after taking a{t) quadratic in t, for instance, a{t) = aot"^, with a^ = 5/27r, 60 = 3 and 
cosx = 4/5. This choice immediately satisfies the loop condition a{T = 27r) = IOtt and 
/3(27r) = 27r, and leads to J^"" cos a{t)dt = 7rC{2V5)/V5 = 0.700896 ^ 0, where C{x) is 
the Fresnel cosine integral. 

Up to now, we have seen that the loop conditions (|3.8| ) guarantee the cyclic time 
evolution of any initial state (|3.5| ) (or equivalently ( |3.7|) ). It should be clear now that 
without these restrictions on a{t) and (3{t), an arbitrary state not necessarily will be 
cyclic. An interesting question arises: are the loop conditions ( |3.8| ) the only way to 
ensure cyclic evolutions of the involved states? In order to get an answer let us consider 
two special initial states of the spin-1/2 system. Let us make in (3.7) 00 = and ^o = X? 
and denote the resulting vector by n+(0); now let us make 0o = ^i" and 9q = n — x, 
and denote the resulting vector by n^(0). The corresponding evolution ( ^.6| ) leads to 



n^(t), where n^(0) = ±e^- Hence, if the spin state points initially along ±e^, the 
rotataion around that vector has no effect on it, and we have nJ:(t) = sinxcos/3(t)i + 
sinxsin/?(t)j + cos^k, n^(t) = ~^x^'t)- At times T„, such that /5(T„) = 2mi, these 
states turn back to n^(0). The subtended solid angles are Afi"^ = 2n7r(l =1= cosx), and 
so the geometric phases become: 

7± = -n7r(l^cosx)- (3.12) 

Let us notice that n represents here the number of effective turns that nj(t) performs 
around the z-axis, and that (|3.12|) is a general result vahd for any x, with a{t) and 



(3{t) arbitrary. In particular, when x = ^/2 (orthogonal rotations) we get 7 = —mr; 
this phase is zero (modulo 27r) when n is even, while it is vr (modulo 2tt) when n is 
odd. On the other hand, if the number of effective turns is zero {x arbitrary) we get 
7^ = 0. This is indeed the trivial case discussed by Zhang et.al. with a{t) = aot, 
[3{t) = /3osm{ujt) and x = 7r/2 §. 



In conclusion, we have shown that the fields ( p.4|) , derived by using inverse tech- 
niques in p, produce exactly solvable spin-1/2 Hamiltonians which, at first sight, should 
be solved by approximate or numerical methods. We have shown also that the time 
evolution induced by these Hamiltonians can produce evolution loops, which simplify 
the calculation of the integrals involved in the formula for the corresponding geometric 
phases. 
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